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Semi-classical Dirac vacuum 
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We study vacuum polarisation effects of a Dirac field coupled to an 
external scalar field and derive a semi-classical expansion of the regu¬ 
larised vacuum energy. The leading order of this expansion is given by 
a classical formula due to Chin, Lee-Wick and Walecka, for which our 
result provides the first rigorous proof. We then discuss applications to 
the non-relativistic large-coupling limit of an interacting system, and to 
the stability of homogeneous systems. 

1 Introduction 

The three-dimensional, massive Dirac operator in a scalar field (j) is given by 
3 Q 

D^ = oij-^ + /3(1 + b) , 

j=i i 

where aj,f3 G are a representation of the anti-commuting Dirac-matrices with 
Oj = 1, 4^ = 1 and we have chosen units such that h = c = m = 1. Under 
appropriate assumptions on (/> it is a self-adjoint operator with domain C^) C 

L^(M^,C^). It arises, for instance, in nuclear physics, where the interaction between 
nucleons is solely modelled by various meson fields, including the scalar fi-meson cj). 
According to the picture of the Dirac-sea, a particle, such as a nucleon, should be 
described by a function ■0 in the positive spectral subspace for Dip in combination 
with the filled ‘sea’ of negative energy states. The state of this particle is then 
described by its one-particle density matrix 

7 := X(-oofi]{Dip) + |0)(0| , 
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where xi denotes the characteristic function of I. In this picture, the spectral 
projection X(-oo,o](-^</>) describes a state with no particles, the vacuum. The total 
energy is then (formally) given by 

Tr(L»07) = -^ Tr \D^\ + (V’, D^, . 

The first term is the formal expression for the energy of the vacuum since the 
spectrum of is symmetric with respect to zero. Of course, — (1/2) Tr is 
not finite and must be regularised. 

Similar objects arise if one considers instead of the scalar field /3(/> an electrostatic 
field (j). Using various methods of regularisation, these have been studied rigorously 
by many authors, especially from a variational point of view o n n osi m 151 
[nmu]. The charge density of the vacuum was derived in a perturbative regime by 
Hainzl jl,‘l| . 

Of the several possible regularisation procedures we will follow [51 EH ED] and 
define the regularised vacuum energy by 

Tvac(</-)=TV(-i|D^|-i?(</-)) , (1) 

where i?((/>) is the fourth order Taylor expansion of —at (/ = 0. It will be 
proved in Proposition |2.l| below that — R{4>) is trace-class under the sole 

assumption that </> £ (M^, M). This is a reasonable assumption as it it is equivalent 

to the finiteness of the energy of the cr-field 


If we consider the total energy, the sum of the energies of the u-field and the Dirac 
particles, the stationary points solve a system of coupled equations. We discuss 
applications of our results, obtained for fixed (p, to this problem in Sections 4.1 
The dynamics of such coupled systems were studied by Sabin 


and 4.2 


with a 

scalar field solving a non-linear wave equation, and in HSIES] for an electric field. 

The formula 0 is very complicated to evaluate, even numerically, for a given (p. 
In applications it is therefore important to obtain simple approximations of Tvac in 
different regimes. For the study of nuclear matter, Walecka [571 ED], Lee-Wick EH 
p. 2306] and Chin (HI Eq. (3.53)] calculated the vacuum energy per unit volume for 
a constant p £ (~1, +oo) and obtained the simple formula 


V(((>) = -^(1 + '/>)'^log(l + p)- P{p ), 


( 2 ) 


where 

P{P) 


1 

4tt‘^ 


7 o 13,3 

7 + 2 ^+^^ + 
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is the fourth-order Taylor polynomial of the first term. In this paper, we give the 
first rigorous derivation of the effective vacuum energy Q , in a semi-classical regime 
where the field 4> varies slowly. Our main result is the following 

Theorem 1.1 (Derivation of the effective vacuum energy). Let 4> G be 

sueh that 


• {(Pl < 1 ; 

• (1 -h |x|2)(i+l“l)/2a“(/> E for every multiindex a E N^. 

Then we have 

lini e^£ivac (</>(£•))= / V(((>(x))dx, (3) 

7r3 

where V is given by (©• 

The convergence (|^ holds under weaker assumptions on cp, but in Section we 
actually prove a more precise result. Namely, we show the existence of a function 
Vs which approximates the true vacuum energy to any order, and which coincides 
with V to first order. That is, 

e^£ivac(</>(£•)) = [ Vs{(p){x)dx + 0{e'^), Vn > 1 


where = V -|- 0{s‘^). The function Ve is a series in e, each term depending on 
a finite number of derivatives of (p. We thereby justify as well the higher-order 
expansions found in later works [21 El El ES]. Explicit formulas of the type (|^ can 
be obtained for all odd dimensions. We will focus on the most interesting case of 
three-dimensional space, the corresponding proofs for other dimensions can easily 
be obtained by adapting our arguments line-by-line. 

As we said, the regularisation used in 0 is not the only possible choice. Follow¬ 
ing one can also consider the Pauli-Villars energy 

which corresponds to introducing new masses mj and coefficients Cj such that cq = 
mo = 1 and 

4 

^Cj(mj)^ = 0, V/c = 0,..., 3. 

j=0 

Our method can be generalized to this case, leading to the convergence 

lime^Tpv(<(>(£•)) = 
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This model still depends on the regularisation masses. In the limit where rrij —)• oo 
for all j = 1, ..,3, one exactly recovers V((/>). 

The advantage of the Pauli-Villars scheme is that it preserves gauge invariance. 
Since the cr-model studied in this paper has no such invariance, the simpler regu¬ 
larisation procedure 0 works just as well. Using the Pauli-Villars method will be 
important for the study of the Dirac vacuum in a slowly varying magnetic field in 
the forthcoming work m- In this case, the effective energy obtained in the limit is 
the famous Heisenberg-Euler Lagrangian EHIEIIES]. 

In physical quantities the parameter e may correspond to h/mc, where m is the 
mass of the Dirac particles, in which case the limit e —>■ 0 would be relevant to 
both the semi-classical limit —?■ 0 and the non-relativistic limit c —oo. For an 
interacting system, for which bound states exhibiting a certain scaling behaviour 
may exist, other identifications are possible. As an application, we discuss the 
simultaneous non-relativistic and large-coupling limit of the Dirac-sigma model in 
Section 4.1 where e is y/h/mc. Finally, in Section 4.2 we present some numerical 
results about the stability of a homogeneous system of Dirac particles interacting 
with a scalar field, a question related to nuclear matter. 


2 The regularised vacuum energy 

From now on we always work with the Dirac operator including a semi-classical 
parameter e 


d 

:= + /3{1 + (p) , 

j=l 3 

which is of course equivalent to taking a slowly varying field (jyeix) = by a 

unitary dilation. In this section we introduce the regularised vacuum energy and 
derive its basic properties. We suppose p G which ensures that ficj) is 

infinitesimally Do-bounded and thus that this operator is self-adjoint with domain 
Dom(D0) = C^). We define the regularised vacuum energy by the formula: 


£va,c{P) • — 



3=0 



(4) 


Note that for small s, zero is not in the spectrum of so we can expect that 
\Dgfp\ is differentiable. Our proof that 0 is a correct definition relies on the integral 
formula 


= 11 
vr Jr 


Dl 


Dl+U3‘^ 


do; = — 
27r. 


2 - 


lUJ 


+ 


ICO 


Ds -I- iw Ds- ic3 


dcj, 


(5) 
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where the integral (in the sense of Bochner) exists in the norm topology of L^), 

even if zero is in the spectrum of D^, because in that norm the integrand decays 
like uj~‘^ for w —)■ oo and stays bounded at w = 0. 

Proposition 2.1. The regularised vaeuum energy defined by equation fl is a con¬ 
tinuous function of 4> ^ It equals 


Svsicifi) = - — Tr 

ZTT 




Dl + u^ 


and is smooth on the open set U := 




^s,p 


s=0 ^S(l> 


DL + uj^ 





Before proving Proposition |2.1| let us prove two preparatory lemmas on the re¬ 
solvent of D^. Throughout, we will denote by I < p < oo, the p-th 

Schatten-ideal in the bounded operators on the Hilbert space J^. The most 

important cases are the trace class {p = 1) and the compact operators {p = oo). 


Lemma 2.2. For every a; G M and 3 < q < oo the map 4> fif’iDo + ^ is 

continuous from H^{R^,R) to ©q(L^(]R^, C"^)). 


Proof. Using the Kato-Seiler-Simon (KSS)-inequality j28[ Theorem 4.1] 


\\f{x)9{-^'^x)\\6fiL'2{R3^) < (2vr) ^^'^||/||l<?(r3)||5||l9(m3) , 


( 6 ) 


and the Sobolev embedding H^{R^) we have for 3 < q < 6: 


with some constant C > 0. Since Sr &q continuously for r > q this proves the 
claim. □ 


Lemma 2.3. The set U := {fi £ H^{R^,R) : 0 ^ a{D^)} is an open subset of 
H^{R^,R) and for every co G M the map i-)- (H^-l-ia;)”^ from U to .if C^)) 

has infinitely many continuous derivatives. 

Proof. If (j) £ U and f £ H^{R^) has sufficiently small norm, then has 

operator norm less than one. Hence -|- fif, is invertible by a Neumann series and 
U is open. Then, the resolvent formula 

{D^ + ia;)-i - {D^ ± = -{D^ + iw)-± iw)“^ (8) 

holds and implies continuity. Letting f, = sf this allows us to calculate the derivative 

(Diff((D 0 ± iw)“^), C) = lim I [{D^ -£ fisf ± iw)"^ - {D^ ± iw)”^) 

= -(D0 ± ± iw)“^, 
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which is obviously also continuous. Using Q again shows that convergence to the 
derivative holds also in the operator norm, and the fact that the derivative is a 
product of differentiable functions gives smoothness. □ 


Proof of Proposition 2^. Consider \Ds(j,\ given by equation (|^. By Lemma 2.3 the 


integrand is a smooth function of s in a neighbourhood of s = 0. One easily checks 
that its derivatives at s = 0 are elements of ((M,dtu),L^)), so we can 
exchange integration and differentiation w.r.t. s, which establishes equality of the 
formula for £ivac given in the Proposition with Q, in the sense that the expressions 
inside the trace are equal. In order to show that this expression is an operator of 
trace-class we first calculate 


{Ds^±iuj) ^ = {-ly jl {Do ± iuj) ^ {/34>{Do±iuj) . 

s=0 

Then, using (§ and the resolvent formula Q , we can write the sum of the integrands 
as ^{I{uj, (f>) + I{—(jJ, 4>)) with 


ds-? 


^ 1 d-^' 

I{uj,f)) : = -iu}{D^ + 

j=0 ^ 


s=0 


{Dstj) -|- 


\UJ] 




= \u{D^ + \uj) ^{l3(f{Do + iu) 

Using the estimate 0 with g = 5 we find 

\\I{oj,(j))\\ei < \\^^{D4, + iu})~^\\e^\\\{P4>{Do + iuj)~^\^\ 
<C{l+uP‘) • 


(9) 


6i 


( 10 ) 


This shows that I{uj,(j)) £ L^(M, ©i) for all (f £ so the expression Q is indeed 
well-defined. 

Continuity in cf is an immediate consequence of Lemma |2.2| and the Dominated 
Convergence Theorem. Since (Diff(,0(/)(Do-|-ia;)“^)), .^) = /3^(Do-|-iw)“^) the deriva¬ 
tives of I{uj,4>) also satisfy an estimate of the form (10), by the Holder inequality 
for traces. Hence smoothness of Tvac also follows from dominated convergence. □ 


3 The semi-classical expansion of the vacuum energy 

In this section we derive the formula ([^ as the leading term in an expansion of 
Tvac in powers of e. The semi-classical expansion is naturally also an expansion in 
orders of derivatives acting on 4>. As such, it reproduces the ‘derivative expansions’ 
discussed in the physics literature [211510123]. For this reason we make very strong 
smoothness assumptions on the field 4> for its derivation. These can be somewhat 
relaxed if one is only interested in a specific order of the expansion. 


6 







We will also need to assume that |(/>| < 1. The lower bound i;/) > — 1 is certainly 
necessary due to the logarithm in formula (§, and the upper bound is needed to 
obtain the following series expansion: 

f Tr = /tt — {/3(j){Do + duj . (11) 

J J j=5 

To prove this formula, we use that \\/3(I){Dq +iw)~^||eoo < 1 write 


Tr [I{uj, (j))] = iwTr 
= iwTr 


{1 + P(j){Do + iuj) ^{Do + iu}) ^{(34>{Do + iu}) 

{Do + iu})~^{P(j){Do + iw)“^)^(l + l3(f>{Do + iw)“^)“^ 


= iw^(-l)^Tr {Do + ioj) ^{/3cj){Do + ioj) 
j=o 


-i\i+5 


Then we note that 

Tr (/3(^(T>o + = -iujj Tr {Dq + iw)"^ (/3<('(T>o + iw)“^)^ 

and integrate by parts, which is justified by the KSS-inequality ([^ and gives ( [IT] ). 
Now ~ log(l + x) — r{x), with a polynomial of order four r{x). As 

I{—u}) = I{oj)* we then find 


1 

Tvac(0) = 


Tr 
+ Tr 


log (1 +/30(T>o + iw) y+r{P(t){Do + \uj) 
log (1 + {Do - \ijj)~^l34)) + r{{Do - ia;)“^/3(()) 


( 12 ) 


do;. 


Under sufficiently strong assumptions on (j) we can represent this quantity as an 
asymptotic series in e, using techniques of semi-classical pseudo-differential calculus. 
We will mainly refer to uni for results on this topic. This book treats only scalar¬ 
valued symbols, as opposed to the valued symbols we use, so we will indicate in 
which sense the results there translate to our setting, with only minor modifications 
of the proofs. 

Theorem 3.1. Let cj) G be sueh that for every multiindex a G there exists a 
eonstant Ca, with cq < 1, satisfying 


9l“l 


(j){x) 


< Ca{l + |x|^) 


(l+l«l)/2 


Then ^vac(^) admits an asymptotic expansion at e = 0 of the form 

. n 

e^£^^c{4>) = T e^Tfc((/))(x) dx -g C>(e"+^), 

JK3 


k=0 
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where Vn{4’){x) depends on the derivatives of order at most n of (p at x and 


= -^(1 + (p{x)flog{l + (p{x)) - P{(p{x)) 

Vimx)=0. 


Proof. In order to keep track of the dependence on to we will use a pseudo-differential 
calculus depending on uj as an auxiliary parameter. That is, we define the symbol 
class 5™ to be those a G x x M, satisfying 


da g/S 


dx°‘ dph 


a{x,p,uj] 


£4x4 


< 


m—\a\ 

Cas{} + \xf) ^ [l + \p\^ + 


2 


for all multiindices a,f3 G and appropriate constants Coi,i 3 . An e-admissible 
symbol of order m G Z is then a function e e-)- a{e) in ^“([0, eo),5™) with aj\f := 
le— 0 ® ^ S™'~^, for which the remainder term of the Taylor expansion of order N 
at e = 0 is of order in the Frechet topology of Let A = f34>{x){Do -|- 

ia;)“^. Because {Dq + ia;)“^ is the quantisation of {a ■ p + (3 + icu)”^ there is an 
e-admissible symbol a(e) of order —1 such that A equals the e-Weyl-quantisation 
Op^ (a), i.e. 

{Af){x) = + y)/2,p,u},e)f{y)dydp, 


for every Schwartz-function / (as follows from inspection of m Proposition 7.7]). 
The leading term is 


ao = pp{x){a ■ p +(3+ iuj) ^g5 


Because ||0||oo < cq < 1, A has norm less than cq. Since F{z) = log(l + z) — r{z) is a 
holomorphic function on the open unit disk, F{A) can be computed using Cauchy’s 
formula. Let j := {z : \z\ = (1 -|- co)/2}. Then we have for z G 7 

-1 2 

\\{ao{x,p)-z) ||c4x4 < --, 

1 - Co 

so {ao{x,p) — z)~^ G 5*^, with bounds independent of z. This implies that {A—z)~^ = 
Op^ {b{z)), for an admissible symbol b{z) of order zero, with bo{z) = {ao{x, p) — z)~^ 
(cf. |ini ( 8 . 11 )], note that the correct generalisation of a lower bound on a scalar 
symbol is an upper bound on the inverse). Thus 

F{A) = ^OpY ^F(z) 6 (z)dz) , 
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so F(A) itself is the quantisation of an admissible symbol / of order zero. Its 
expansion is given by 


/n = ^ y F{z)bn{z) dz , 

in particular /o = F(ao) (the formulas for the higher order terms will differ from 
those for scalar-valued symbols). We now use this to evaluate the trace and the 
integral over oj, order-by-order in n. For n = 0 we have 

/o = F(ao)=a5^^a^oG‘5-^ 
i=o ^ ^ 

since the sum is absolutely convergent. This implies that /o and its derivatives are 
elements of so (cf. [IHl Theorem 9.4]) Op^ (/o) G Si and 

Tr[Op^ (/o)] = ttt'ws / ^^'c4fo{x,p)]dpdx . 

(27re)-^ Jr6 

The leading order in the expansion of V is thus given by 

Vo{(j}) = - '3 [ f Trc4 [F(ao(x,p))] + F{ao{x,p)*)] dpduj 

47r{27ry^ Jrs 

= - ^^^27r)3 / 

Now, 


(1 -|- ao)(l -|- Oq) — 1 -|- uq -|- Og -|- ogag 

= 1 -h {2(j){x) + (j){x)‘^){p"^ -h 1 -h u }^)~^, 

so the C'^-trace is trivial and expanding into power series and setting ^ = (p, uj) G 
we obtain: 

= f f:^^^(e' + l)-^■(20(x) + </>2(x))^■-r'((/)(x),Od^ 

JK4 ^ J 


The new remainder term r' is a polynomial of fourth order in (^{x ), and it is integrable 
in since the first term and the total integrand both are. The ^-integral of the first 
term evaluates to 


Jr* 


0 --l)(j- 2 )’ 
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and the sum then yields 


OO 


E 

J=3 




{2(j){x) + (i)‘^{x)y 


= (1 + 2(j){x) + 0 ^(x)) log(l + 2cj){x) + (/)^(x)) - r''{(j){x)) 
= 2(1 + 0 (x))^log(l + (j){x)) - r"((/>(x)) , 


with another fourth-order polynomial r". With this we obtain 

Vo(<(') = + (t>{x)y\og{l + (t>{x)) + r"{(l){x)) + • 

'-V-' 

= :-Pi<f>(x)) 

By construction, P is a fourth order polynomial and since the first four derivatives 
of the integrand in equation (13) w.r.t. cj) vanish, P{t) is the Taylor polynomial of 
-4^(1 + t)'^log(l + t). 

We have thus found the claimed leading order term for Tvac(0) and it remains to 
estimate the error. First, note that / — J2n=o^^fn ^ and that this implies a 
bound in L^(M,©i) on the quantisation, because 

Opf ((1 + x2)-5/2(i + p 2 ^ a;2)-5/2) g l\R, © 1 ), 
and for every cj G 5“^ 


OpJ'' ((1 + Opf (a) 


<C. 


by the Calderon-Vaillancourt Theorem m Theorem 7.11]. This proves that the 
expansion of / gives an expansion of £va,c{4’) of the claimed form, with 


Vn{^)ix) := [ [ TTc4[fnix,p,u})]dpduj, 

JR JR^ 

for the terms of order at least five. It thus remains to estimate fn for 1 < n < 4. For 
brevity, we will only perform the proof for fi , as the arguments for higher order terms 
are essentially the same. The explicit form of 61 , and thereby /i, is obtained from 
the formula for the Weyl-Moyal product (see cni Theorem 7.3]), which represents 
the product of operators on the space of admissible symbols. Denoting this product 
by o and by {•, •} the Poisson-bracket we have 


bi{z) = -bo{z)aibo{z) - {bo{z) o (ag - z))ibo(z) 

= -(ao - z)~^ai(ao - z)~^ - |{(ao - 2 ;)“\ao}(ao - z)~^ 
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Now for every a G S'^ we have 

[a, (ao - z)~^] = -{ao - z)“^[cr, ao](ao - z)~^ G 5™“^ , 
so the first term of 6 i, for example, can be written as 

bo{z)aibo{z) = 0160(2;)“^ + [oq, ai]5o(2;)“^ + [oq, [oq, ai]]bo{z)~‘^ + S~^ 
On the other hand, we have 


1 

27ri 



.(- 1 ) 


1 


i-i 


n 

\e-i 


{£-l)ldz^ 
d^-i 


ct(«o - z) 


-1 


dz 




and since 


(^f) (oo) G 5 for £ < 5 this implies that 


1 

27ri 





ai(ao 


z)-^dz G 5“^ 


The same argument, applied to the second term of 61 , gives /i G S~^, which proves 
the desired bound. 

To complete the proof, we still need to check that Vi = 0. We have 


TV, 


'C4 


— y F( 2 ;)(ao - 2;) ^01(00-2;) ^dz 


Trc4[aiF'(oo)] 


and 


ai{x,p) = {/34>{x)o{ap + /3 + iuj) = -^{Pcl),{ap +(3+ iuj) ^} = -^yxypao- 

Thus, for every x with (l){x) 7 ^ 0, 

TVc4[aiF'(ao)](x,p) = -5VpTrc4[(Vlog(/)(x))F(ao(x,p))] , 

the dp-integral of which vanishes. For the second term of 61 , the C'^-trace equals 
zero, because V^ao commutes with ao and thus 

Trc 4 ((Va;(ao - 2 :)“^)(Vpao)(ao - 2;)“^) = 

= - Trc 4 ((ao - 2;)“^(Vxao)(ao - 2;)“^(Vpao)(ao - 2;)“^) 

= Trc 4 ((ao - 2;)“^(Va,ao)(Vp(ao - ^)“^)) 

= Trc 4 ((Vp(ao - 2;)“^)(Vxao)(ao - 2;)“^)) . 

This concludes the proof of Theorem |3.1| □ 
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4 Applications 


4.1 An interacting Dirac Klein-Gordon model 

In this section we discuss the implications of the results of Section for a simple 
model of a particle, described by a Dirac equation, interacting with a classical field, 
described by a Klein-Gordon equation. Such models arise in the mean-field approx¬ 
imation of models in nuclear physics. They are often treated with the additional 
‘no sea’ approximation, which amounts to neglecting the contribution of negative 
energy states. In such models, there are usually additional vector-fields to consider. 
However, due to gauge-invariance of the Dirac equation, the leading order of the 
vacuum energy is solely due to the scalar field. With the ‘no sea’ approximation, 
the static model amounts to finding critical points of the energy 


To (V’,0-) 

:= [ 'ip{x) (—ica ■ V +/3{c^ + ga{x))'] + lc^\Va{x)\‘^ + lM‘^c^\a{x)\'^dx , 

under the constraint that 'll; G be a normalised function in the positive 

spectral subspace of —ica • V-|- /3(c^ -|- g4>{x)), where we have introduced the physical 
constants, c (speed of light), M (mass of the field), and g (particle-field coupling 
constant). If we rescale lengths by e = ^/c, choose large couplings of the form 
g = goe~'^^'^ and set (f){x) = goa{£x) this becomes 

To (V', 4>) = J^ tpix) (^-iea - V + /3{1 + £^4>{x))^ '4>{x) 

+ 55o (e^|V(/)(x)|^ -h M‘^\(p{x)\‘^'^ dx . 

Writing ?/) = (y?,y) with ip,x ^ i?^(M^,C^), the Euler-Lagrange equations with a 
Lagrange multiplier 1 — e'^g, are given by 

{ —iea ■ Vx = 

-iecT ■Vg= {2- g + £^4>ip)x 
-e^^c^ + M^=-gl{g^-X^). 

For £ small enough, these equations are known [221 Theorem 3] to have a solution 
{ips,Xei4’e)^ for which —)■ u^/’nls fo T1^(M^,C^), where u G C is an arbi¬ 
trary vector and V’nls is the unique positive solution to the non-linear Schrodinger 
equation 

- (11)^ IV’lV + l^V’= 0. (14) 

By a scaling argument, there exists a unique g for which ||'0 nls||l 2 = 1, and the 
corresponding solution Xe, (j^e) is a critical point of Tq satisfying the constraints. 
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If we apply the same procedure to the energy £ including the contribution of the 
‘sea’ of negative energy states, defined as in Section we obtain 

£{ii, (j)) = Soii/j, 4>) + e“^£ivac(e^0) • 


Remember that we have shown in Proposition 2.1 that £va,c{<P) is smooth in an H 


neighbourhood of (/> = 0, and that its first four derivatives at this point vanish. 
Hence the contribution of will be small, of order for e —)• 0, despite the 

pre-factor of £~^. It is thus easy to show, by the same implicit-function argument 
used in [22], that this energy also has a critical point exhibiting the same 

limit, as e —)■ 0, as without the vacuum contribution. For the term £ivac(£^</>) both 
an expansion around (p = 0 and the semi-classical expansion of Section are appro¬ 
priate. The second is justified by the fact that the solution to the NTS equation (141 
is smooth and exponentially decreasing (see j29|b so the hypothesis of Theorem 3.1 
are fulfilled in this limit. 


4.2 The stability of homogeneous systems 

In this section we discuss the effects of vacuum polarisation on the stability of 
homogeneous systems. For uniform nuclear matter, stability was already studied in 
the early works laiziis] and many later contributions, but with an emphasis on the 
dominant effect of the vector mesons and neglecting vacuum polarisation. 

In this section we present formal derivations and numerics, we will indicate how 
some of these can be made rigorous. We take the state of the system to be of uniform 
density and chemical potential /x > 0, 

T’m X(—oo,/i] • 

The energy of the positive energy states is then formally given by 

where M is the mass of the field and g the coupling constant. The total energy is 
the sum £o[g, </>) +£va.c{(p)- For a constant field, the energy per unit volume becomes 


252!'^!^+ (27r)3 


+ {I + (pydp -h v{4 >), 


(15) 




with V{(p) given by equation @- 

We define (^vac £ to be the local minimum of (15) in this interval and 

(po G (“1)0] to be that of T — V (see Figure [^. Note that (()vac is only a local 
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minimum, as V is unbounded from below for positive values of (j). We will now study 
the response of the field (p^ac to a local perturbation r] of this system by calculating 
the response function T, which is formally given by 

<?(A^,</>vac + <5?/) = 77 ^ / |r}(p)pr(p,/i)dp, 

5=0 ( 2 ^) -'R3 

and compare this to the corresponding results for (pQ. This expression is only formal 
since (pvac + 6r] ^ and £vac{(pvac + St]) is not defined. In order to give a proper 
definition, one should consider a system in a finite box and take the thermodynamic 
limit of the difference £{n, (pvac + ^v) (pvac), with r] G H^. We will not perform 

this construction here, the function T will give a finite expression on the right hand 
side for sufficiently regular r]. 

The response of the vacuum can be obtained by derivation of the formula (cf. ( |12[ ) ) 

£vac{(p+Sv) = - 7 !- / Tt (log {l + 5Pr]{Do +P(p+iuj)~^)-r{(p+6r],u:)) duj. (16) 
27r ^ ^ 

Calculations similar to those of m Lemma 4.2], the details of which are given in 
the Appendix, give the result 

T^acip) =(27r)^V"(0vac) + ( u{l - u) log i du 

40 \{l + (p^s-cY + u{l - u)p^ J 

r, 2 (pvacu{l — u) / 2 o22|ft32 

+ 27rp / 2 - 6tt + 2up - 8u p + 6u^p 

Jo (1 + u(l - ujp^Y V 

+ 0vac(l + 2ttp^ — + lOu^p"^ — 5u^p‘^)'^ du . 

The complicated form of this expression is largely due to the greater number of 
regularising terms as compared to m- The response of the positive energy states 
is given by 

- 27rir+(/i,p) 

= f f Trc4 [/3{aq + j3{l + c^ivac) + z)~^f3{a{q - p) + /3(1 + (pvac) + z)~^] dqdz , 

Cm R3 

where is a path in the resolvent set of the matrices ap + /3(l + (pvac) enclosing the 
interval (0, p]. The total response function is then given by 

T{p,p) = ^ + ^ + + ^vac(p) ■ 


d^ 

d^ 
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Figure 1: 0vac 
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and kp = 
(respectively). 


(solid line) 
line) as func- 
V/r^-(l + <(>0)^ 

(/)vac)^ 



Figure 2: Response functions T{n,p) relative to (dashed line, without vacuum 
corrections) and (p^ac (solid line, including these) as functions of kp, at fixed p = 0.01 
(left) and p = 1 (right). 


Numerical evaluation of these quantities, introducing the Fermi momentum kp := 
— (1 + (pvac)'^ and choosing the parameters g = 10 and M = ^ (this choice is 
suggested by [Ml Table 3]), gives the plots shown in Figure]^ 

We observe that the vacuum polarisation has a twofold stabilising effect: First, 
it stops (p from reaching the point —1 where the gap in the spectrum of D(j) closes, 
and second it renders the minimizer cp more stable to local perturbations compared 
to 00, which can be unstable at high Fermi momentum. It is an interesting open 
problem to make these observations rigorous. 
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Appendix: Calculation of the vacuum response function 


We start by taking the derivative of the integrand in (161, writing cp = cpy 


d^ 

d(i2 


( 5=0 


Tr log (1 -h 6(3r]{Do + j3(p + iw) ^) - r{(p -h 6 r], u) 
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with 


^ ('—lii+i 
r{(j) + 5t], uj) = J2 - -^- 


j=i 


+ dr]) {Do + iu! 


,-i 


We have 

d52 

and 

d52 


Tr 


( 5=0 


(^log {I + S/3ri{Do 


lUJ] 


\-i' 


= -Tr 


(^r/CDo +/30 + iw) 


Tr (r((^ + 5r], oj)) 


( 5=0 


= Tr 


- (/3ry(T» + + 2(/)/3(D + {^r]{D + iw)-^) 


4, ,5-i\2 


- (j)^[{I3{D + \uj) ^l3ri{D + \uj) + 2(/3(T> + iw) ^)^(/3r?(T) + iw) 


- 1 \ 2 / 


\-l\2 


which is only formal as the terms are not individually trace-class, but can be made 
rigorous by inserting ±f{5ri, u) which is the remainder for a Dirac operator Dq with 
mass 1 + 4>. 

Now let f,g^ x C^) be sufficiently regular and decay at infinity. Then 

the operator r]g{Do + \uj)gf{Do + icj) has a continuous integral kernel and 


Tr {gg{DQ + iuj)gf{Do + iw)) 

= ^^;^yi??(p)I^Trc4 g{aq + l3 + iuj)f{a{q-p) + P + iu}) dpdk . 


(17) 


This gives us (writing f{q) = {aq + P + iw) and g{q) = {ap -|- /3(1 + ({))+ iw) 
Tvac(p) = ^ / Ti-c^ fdg{q)Pg{q -p)- I3f{q)f3f{q - p) + 2(f){f3f (q))^(3f {q - p) 


-r 


\ldf{q - P))^ + ‘2{/3f{q))^/3f{q - p)) dgdw . 


Using that g{q) = {aq + ,0(1 + ({>)— iu})/{q‘^ -|- (1 -|- + oj"^) the spinor trace of the 

first term can be evaluated, yielding 


Trc4 [Pg{q)Pg{q - p)] 


4(1 -|- 0)2 — 4a;2 -|- Aqp — Aq^ 

(g2 -h (1 -h 0)2 0j‘^){q - p)2 -h (1 -h 0)2 -h a;2) ' 


(18) 


The first term involving / gives the same result with ‘mass’ one instead of 1 -|- 0, 
and for the remaining ones we have 


Trc4 [/3/(g)'+W(g-p)'+1 


Oj,k{^,q,p) 

(g2 + 1 + w2)l+i((g _ p)2 + 1 + ^2)l+fc ’ 
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with the polynomials 

9i^q =4 + 8pq - 12q^ - 120;^ 

^2,0 =4 + I2pq — 24q^ — 4(pg)g^ + — 24a;^ + Apquj'^ + Sq^u? + 4a;^ 

01,1 =4 — 4p^ + 24pg — 24g^ + 4p^q^ — 8{pq)q^ + 4g^ 

- 24w^ + 4p^a;^ - Spqu}"^ + Sq^u? + 4uj^ . 


Using the substitution 

CO OO ^ 

^i+fj^i+k = J y d^dC = J y si+J+V(l - n) dsdn 


1 OO 


0 0 


0 0 


leads us to the Gaussian integrals 

[ do;, [ d^. ^ 

for Lo and the three components of q. Evaluating these gives the terms 

with m = l,orm = l + (/) for the evaluation of (181, and 

©0,0 = - f + 4m^ + 4ri(l - u)p‘^ 

01.0 = — ^ + 4 — 4p^ + 16p^n — 12p^u^ 

©2,0 =4 - 48 + 8p2 - 24p\ + 24p\^) 

+ 4 — 4p^ + 24p^u — 24p^u^ + 4p^M^ — 8p‘^u^ + 4p^u^ 

011 =4 + I2p^ - 36p^u + 24p^u^ - 48 ) 

— 12p^ + 36p^u — 4p'^u — 24p^v? + 12p^tt^ — 12^^11^ + 4p^u^. 

After gathering all of these terms (with the correct pre-factors) we split T^adp) = 
Evac(O) +Tsing(p) +ri.eg(p), where Ts ing contains all those summands whose integrand 
vanishes at p = 0 and diverges at s = 0. We then have 

rvac(O) = i27rfv"{cP) , 

since p = 0 corresponds to a perturbation by a constant field. Calculating the ds 
integrals for the other terms yields 


rl poo 

Tsi„g(p) = 271 / / u(l - n)p2e-*“(i-©^’ 

Jo Jo s 

„ n ,, / l+tt(l— u)p^ 

= 27 rp^ J rt(l — u) log ( 


_ -s 


(1 + + m( 1 — M)p2 


• dsdtt 
du , 
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and 


Tregip) = 27rp^ f , 2)2 “ 6u + 2up^ - 8u^p^ + 6wV 

Jo [1 + u{l - u)p^)^ 

+ </)(! + 2 up^ — 7u^p^ + lOu^p^ — dit. 

References 

[1] V. Bach, J. M. Barbaroux, B. HelfFer, and H. Siedentop. On the stability of the 
relativistic electron-positron field. Commun. Math. Phys., 201:445-460, 1999. 

[2] P. Blunden. Evaluation of Dirac sea effects in a finite system. Phys. Rev. C, 
41(4):1851, 1990. 

[3] P. Chaix and D. Iracane. From quantum electrodynamics to mean field theory: 
I. The Bogoliubov-Dirac-Fock formalism. J. Phys. B, 22:3791-3814, 1989. 

[4] P. Chaix, D. Iracane, and P.-L. Lions. From quantum electrodynamics to mean 
field theory: II. Variational stability of the vacuum of quantum electrodynamics 
in the mean-field approximation. J. Phys. B, 22:3815-3828, 1989. 

[5] L.-H. Chan. Derivative expansion for the one-loop effective actions with internal 
symmetry. Phys. Rev. Lett., 57(10):1199, 1986. 

[6] O. Cheyette. Derivative expansion of the effective action. Phys. Rev. Lett., 
55(22):2394, 1985. 

[7] S. A. Chin. Relativistic many-body studies of high density matter. Phys. Lett. 
B, 62:263-267, June 1976. 

[8] S. A. Chin. A relativistic many-body theory of high density matter. Ann. of 
Phys., 108:301-306, 1977. 

[9] S. A. Chin and J. D. Walecka. An equation of state for nuclear and higher- 
density matter based on relativistic mean-field theory. Phys. Lett. B, 52:24-28, 
Sept. 1974. 

[10] M. Dimassi and J. Sjostrand. Speetral asymptoties in the semi-elassieal limit, 
volume 268 of London Mathematieal Soeiety Leeture Note Series. Cambridge 
university press, 1999. 

[11] P. Gravejat, C. Hainzl, M. Lewin, and E. Sere. Construction of the Pauli- 
Villars-regulated Dirac vacuum in electromagnetic fields. Areh. Rat. Meeh. 
Anal, 208(2):603-665, May 2013. 


18 



[12] P. Gravejat, M. Lewin, and E. Sere. Derivation of the Heisenberg-Euler La- 
grangian. in preparation, 2015. 

[13] C. Hainzl. On the vacuum polarization density caused by an external field. 
Ann. Henri Poineare, 5(6):1137”1157, 2004. 

[14] C. Hainzl, M. Lewin, and E. Sere. Existence of a stable polarized vacuum in the 
Bogoliubov-Dirac-Eock approximation. Comniun. Math. Phys., 257(3):515-562, 
2005. 

[15] C. Hainzl, M. Lewin, and E. Sere. Self-consistent solution for the polarized 
vacuum in a no-photon QED model. J. Phys. A, 38(20):4483-4499, 2005. 

[16] C. Hainzl, M. Lewin, and E. Sere. Existence of atoms and molecules in the 
mean-field approximation of no-photon quantum electrodynamics. Areh. Ra¬ 
tion. Meeh. Anal., 192(3):453-499, 2009. 

[17] C. Hainzl, M. Lewin, and J. P. Solovej. The mean-field approximation in 
quantum electrodynamics: the no-photon case. Comniun. Pure Appl. Math., 
60(4):546-596, 2007. 

[18] C. Hainzl, M. Lewin, and C. Sparber. Existence of global-in-time solutions to 
a generalized Dirac-Eock type evolution equation. Lett. Math. Phys., 72(2):99- 
113, 2005. 

[19] C. Hainzl and H. Siedentop. Non-perturbative mass and charge renormalization 
in relativistic no-photon quantum electrodynamics. Commun. Math. Phys., 
243(2):241~260, 2003. 

[20] W. Heisenberg and H. Euler. Eolgerungen aus der Diracschen Theorie des 
Positrons. Z. Phys., 98:714-732, Nov. 1936. 

[21] T. D. Lee and G. G. Wick. Vacuum stability and vacuum excitation in a spin-0 
field theory. Phys. Rev. D, 9:2291-2316, Apr 1974. 

[22] M. Lewin and S. Rota Nodari. Uniqueness and non-degeneracy for a nuclear 
nonlinear Schrodinger equation. Nonlinear Differ. Equ. Appl, in press, 2014. 

[23] R. J. Perry. Effects of the Dirac sea on finite nuclei. Phys. Lett. B, 182(3):269- 
273, 1986. 

[24] P. G. Reinhard. The relativistic mean-field description of nuclei and nuclear 
dynamics. Rep. Prog. Phys., 52(4):439-514, 1989. 

[25] J. Sabin. Global well-posedness for a nonlinear wave equation coupled to the 
Dirac sea. Appl. Math. Res. Express (AMRX), 2014(2):312~331, 2014. 


19 



[26] J. Schwinger. On gauge invariance and vacuum polarization. Phys. Rev., 
82:664-679, 1951. 

[27] B. Serot and J. Walecka. The relativistie nuelear manybody problem, volume 16 
of Advanees in Nuelear Physies, page 212. Plenum, New York, 1986. 

[28] B. Simon. Traee ideals and their applieations, volume 35 of London Mathe- 
matieal Soeiety Leeture Note Series. Cambridge University Press, Cambridge, 
1979. 

[29] W. A. Strauss. Existence of solitary waves in higher dimensions. Commun. 
Math. Phys., 55:149-162, June 1977. 

[30] J. D. Walecka. Theoretieal Nuelear and Subnuelear Physies. Imperial College 
Press and World Scientific Publishing Co. Pte. Ltd., second edition, 2004. 

[31] V. Weisskopf. Uber die Elektrodynamik des Vakuums auf Grund der Quan- 
tentheorie des Elektrons. Math.-Fys. Medd., Danske Vid. Selsk., 16(6):l-39, 
1936. 


20 



